Low-field octupoles and high-field quadrupoles in URu 2 Si2 
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The recent experimental finding of large-amplitude antiferromagnetism induced by uniaxial strain 
shows that the "hidden" low-field order of URu2Si2 breaks time reversal invariance. We propose a 
new crystal field model which supports T% octupolar order in the low-field phase, and quadrupolar 
order in a disjoint high-field phase. The temperature dependence of the linear and third order 
magnetic susceptibility is in good agreement with the observed behavior. 



\_phase transition of 
Q. Though URu 2 Si 2 



The nature of the T = 17.5K 
URU2S12 is a long-standing puzzle 
was long considered as a "light" heavy fermion system, 
implying that the /-states should be included in the 
Fermi volume, many aspects of the normal state behav- 
ior are well described in terms of a localized /-electron 
model. Specific heat measurements [2j show that an elec- 
tronic entropy of 0(\n 2) is released by the time the tem- 
perature reaches 30K, and a sizeable fraction of it is as- 
sociated with the A-anomaly at 17. 5K. Thus the phase 
transition should be associated with the full-scale order- 
ing of a localized degree of freedom per site, but the na- 
ture of the order parameter remains hidden. It is obvi- 
ously not the tiny (M z ~ 0.03/iR, z being the tetragonal 
fourfold axis [001]) antiferromagnetic moment observed 
by neutron scattering [3j. In fact, the association of mi- 
cromagnetism with the 17. 5K transition is dubious, since 
it depends on sample quality, while the thermodynamic 
transition itself is a robust sharp feature. 

If one assumes that the observations are made on 
single-phase specimens, then the weak antiferromag- 
netism should be described as a secondary phenomenon 
driven by the primary ordering of the hidden order pa- 
rameter. The staggered dipole moment and the hidden 
order parameter would possess the same spatial and time 
reversal symmetry. The finding of a first order transition 
to a high-pressure phase with large moments was argued 
to favour this scenario 0. A general symmetry analysis 
listed local octupoles as well as triple-spin correlators . 
We note, however, that an extensive mean field study by 
Santini et al considered the possibility of octupolar order, 
and discarded it in favor of the alternative of quadrupolar 
order Q. Most recently, unconventional density waves 
with alternating plaquette currents were advocated 0. 
The plaquette current can give rise to weak orbital an- 
tiferromagnetism, and is thus in principle well suited to 
describe URu2Si2- However, these works are based on the 
extended s-band Hubbard model; we prefer a description 
emphasizing the orbital character of /-electrons. 

We propose an alternative scenario in which hidden or- 
der and antiferromagnetism are of different symmetries. 
There is experimental evidence that the apparently tiny 
moments belong to a minority phase, and the hidden or- 
der of the majority phase is non-magnetic pi. lllj. 



Assuming that hidden order is not of the same symme- 
try as J z dipoles, there are still two basic options: hidden 
order may, or may not, break time reversal invariance 
Earlier, both possibilities seemed open [lCj, and Santini's 
quadrupolar model gives an example of time reversal in- 
variant hidden order. However, we are going to argue 
that recent observations necessitate to postulate time re- 
versal invariance breaking hidden order. 

According to a crucial recent experiment [Tlj . uniax- 
ial pressure applied in [100] or [110] directions induces a 
relatively large magnetic moment in direction [001]. In 
contrast, stress in the [001] direction does not induce any 
sizeable moment. Since stress is a time reversal invariant 
perturbation, it can induce magnetism only from an un- 
derlying (hidden) order which itself breaks time reversal 
invariance. The directionality of the effect indicates time 
reversal invariance breaking orbital order, most straight- 
forwardly octupolar order. 

Octupolar order as primary order parameter was sug- 
ested for Phase IV of Cei-^LaxBg |l2|, and for NpC>2 



[13J. Knowing such precedents, it is a plausible idea 
to check whether a difficult-to-identify order of an /- 
electron system is octupolar. 

The tetragona l sy mmetry classification of the local or- 
der parameters [lj] in the absence of a magnetic field 
is given in Table |U g and u refer to "even" and "odd" 
under time reversal. Here we neglected wave vector de- 
pendence; it will be specified when we mean staggered 
rather than uniform moments. 



TABLE I: Symmetry classification of the local order parame- 
ters for H = (2?4h notations|ll^. overline means symmetriza- 
tion [3). 
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We argued that the H = order parameter must be 
one of the u operators; it cannot be Ai u or E u because 
that would mean magnetic order; later we mention why 
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it cannot be Ai u ; so it must be one of the octupoles Tf 
(B2u), or 7^1,2 ( B\ u ). Lacking a microscopic analysis of 
the multipolar interactions in URu2Si2, we cannot decide 
between the two, and arbitrarily choose the Tf octupole 
as the zero- field order parameter fl6| . 

Switching on a field H || z, geometrical symmetry is 
lowered to C 4 h- However, the relevant symmetry is not 
purely geometrical. Though taken in itself, reflection in 
the xz plane a ViX is not a symmetry operation (it changes 
the sign of the field), combining it with time reversal T 
gives the symmetry operation Ta v . x . The same holds 
for all vertical mirror planes, and Ci _L z axes, thus the 
full symmetry group consists of eight unitary and eight 
non-unitary symmetry operations |l5| 

Q = Cih + Ta ViX Cih ■ (1) 

We may resort to a simpler description observing that 

Q = C 4 + fa v , x C 4 (2) 

is an important subgroup of Q, and we can base a sym- 
metry classification on it. The multiplication table of Q 
is the same as that of C 4v , and therefore the irreps can be 
given similar labels. It is in this sense that the symmetry 
in the presence of a field H || z can be regarded as C 4v 
(a convention used in ^H). The symmetry classification 
of the local order parameters valid in H || z is given in 
Table ITT1 The results make it explicit that the magnetic 
field mixes dipoles with quadrupolcs, quadrupoles with 
certain octupoles, etc. 

TABLE II: Symmetry classification of the lowest rank local 
order parameters for H || z (notations as for C 4v iJJj) 
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1, J z 
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In a field H || z, there can exist ordered phases with 
four different local symmetries: A2, B\, B2, and E. 
The zero-field i?2«-type Tj 3 octupolar order evolves into 
the -Bi-type T^-0\ mixed octupolar-quadrupolar order 
(Figure [IJ. Experiments tell us that B\ u (and also B\) 
order is alternating (Q = (111)). The gradual suppres- 
sion of octupolar order under field applied in a high- 
symmetry direction is well-known, e.g., from In 
our calculation, the octupolar phase is suppressed at 
H crA re 34.7T (Figure 0. 

There can be other kinds of order, but they cannot 
coexist with B\ because they carry different symmetry 
labels. The phases can be disjoint, separated by non- 
ordered regimes, or when they press against each other, 
the transition must be first order. It is a question of detail 




FIG. 1: The high-field part of the T = phase diagram 
of the multipolar model (H in units of T (Tesla)). Vertical 
axis: (Tf) for the low-field phase, and (O zx ) for the high- 
field phase. The field-induced mixing of the order parameters 
is shown within the shaded areas. The overall appearance of 
the T-H (inset, T in units of K) phase diagram is very similar. 
(The critical temperature of the E phase is scaled up 3-fold). 



whether isolated phases are bounded by critical lines or 
first order boundaries. 

Seeking agreement with high-T and large-field data we 
postulate a crystal-field model in which two levels tend to 
cross at H > H CTi \, and they are connected by matrix el- 
ements of E operators. Consequently, we find a high-field 
E phase where { J x , J y }-type transverse dipolar order is 
mixed with {O xz , Oj,z}-type quadrupolar order (see Fig- 
ure ^and Table ITT|) . The overall appearance of our phase 
diagram closely resembles the results of high-field mea- 
surements |l8j . Some experiments identified additional 
domains in the H—T plane |l9j, but we think that the 
two phases shown in Figure ^ are the most robust part 
of the phase diagram. 

We assume stable 5/ 2 valence, and Hund's rule J = 4 
ground state. Let us seek a plausible level scheme to sup- 
port the postulated ordering phenomena. It is accepted 
that the ground state is a singlet, and it is connected to 
another singlet across a gap of ~ 100K by a matrix ele- 
ment of J z : 3J. In our scheme, \t\) is the ground state, 
and |*2 ) the A2 = 100K excitation. We need the low- 
lying (Ai = 45K) singlet \t 4 ) to allow induced octupole 
order. Finally, as in previous schemes [||, at least two 
more states are needed to fit magnetization data up to 
300K. We found it useful to insert one of the doublets 
(\d±)). Level positions were adjusted to get good overall 
agreement with observations but we did not attempt to 
fine-tune the model. The list of the relevant crystal field 
states is given in Table ITTTl (we use a = 0.98, b = 0.22). 
The field dependence of the levels is shown in Fig. [21 

Our crystal field scheme differs in essential details from 
previous ones [|J, but the quality of fits to the suscepti- 
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TABLE III: Tetragonal crystal field states used in the model 
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FIG. 2: The magnetic field dependence of the single-ion levels. 

bility (Fig|3|), the non- linear susceptibility (Fig0J left), 
and the metamagnetic transition (FigE] right) is not in- 
ferior to what was achieved earlier 20]. However, our 
basic argument in favour of the present scheme is not 
that it recovers standard results, but that it allows the 
derivation of the phase diagram shown in Fig ^ We use 
the mean field decoupled hamiltonian 

Hmf = Ai|i 4 )(* 4 | + A 2 |i 2 )(i 2 | + A 3 K)(d a \ 

«=+,- 

-gfi B HJ z + A oct (Tf) Tf - A quad (O zx ) O zx (3) 

where g = 4/5, and the octupolar mean field coupling 
constant A oc t is meant to include the effective coordina- 
tion number; similarly for the quadrupolar coupling con- 
stant Aq Ua d. We assume alternating octupolar order [23j |. 
and uniform O zx order; the result would be the same if 
the high-field quadrupolar order is also alternating. We 
do not introduce 0\ or {J x , J y } couplings, nevertheless 
(0 2 ) ^ in the B\ phase, and (J x ) ^ in the E phase. 

At H = 0, the only non- vanishing octupolar matrix 
element is C = (ti\Tf\t4) « 8.8. Octupolar order is 
driven by the large C: assuming A oc t = 0.336K we get 
the critical temperature T (H = 0) = 17. 2K for T/-type 
antiferro-octupolar order. Using a similar estimate, we 
find A£P w 0.2K for Np0 2 which orders at 25K 0, thus 
the assumed octupolar coupling strength is not unrea- 
sonable |22^ . 

The octupolar transition shows up as a break in the 
temperature derivative (dxi/dT) of the linear suscepti- 
bility (Fig. |3 right). The sign of the discontinuity of 
slope is related to the fact that the critical tempera- 
ture decreases in magnetic field like To{H) ss To(H = 
0)~a H H 2 dHHl. An Ehrenfest relation HH3 connects 
the discontinuity of (dxi/dT) to that of the non-linear 
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FIG. 3: Linear susceptibility per site (in //b/T) on extended 
temperature scale (left), and in the vicinity of the octupolar 
transition (right). The dashed line gives the single-ion result. 
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FIG. 4: Left: Nonlinear susceptibility X3 a t the octupolar 
transition. Right: the magnetization curve at T = (M in 
units of (ib)- 

susceptibility X3 (Fig-E left). 

Up to the vicinity of the t\-d- level crossing shown 
in Fig. |21 field effects can be understood within the \t\)- 
|t4)— 1*2) subspace. Applying H || z mixes \t 2 ) to \ti). 
This has two effects. First, since 0\ connects |£ 2 ) to ^4), 
the order parameter acquires a mixed T!?-0\ character 
(cf. Table [njl. Second, it reduces the octupolar ma- 
trix element, and thereby also To- With the parameters 
given before, the octupolar transition is fully suppressed 
at -ffcr.i = 34. 7T (see Fig.QJ. There is an accompanying 
change in the slope of the magnetization curve which, 
however, is not noticeable on the scale of Fig. 0] (right). 

A basis-independent description of field effects relies 
on the Landau expansion of the Helmholtz potential A 
which yields the magnetic field as a derived quantity H = 
{dA/dJ) E3. A is a sum of invariants. Two terms which 
are important for the present purpose, are contained in 

l{A 2u ®Bi g ®B2u) = ciJ z (O)T/(Q)02(-Q) 

+c 2 J 2 (Q)T/(-Q)O 2 2 (0). (4) 

ci and c 2 are non-zero even if we consider the lowest 
two levels only. The existence of this invariant can be 
exploited in several ways. In a uniform magnetic field, 
alternating octupolar order induces similarly alter- 
nating quadrupolar order 0\. Alternatively, it follows 
that in the presence of uniform quadrupolar polarization 
02 1 alternating octupolar order gives rise to a magnetic 
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moment J z with the same periodicity. Such a quadrupo- 
lar polarization is created by uniaxial stress applied in the 
[100] direction, which is observed to give rise to alternat- 
ing magnetic moments of O(0.1/iB), clearly different from 
micromagnetism ^lj- Fig. [5] shows the stress- induced 
staggered magnetization for the same set of parameters 
as in previous plots. Sufficiently large stress suppresses 
octupolar order, like a sufficiently strong field does. The 
maximum induced moment is ~ 0.5/ib; the measured 
~ 0.2/ie ^3 m ay belong to the rising part of the curve. 
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FIG. 5: Stress-induced magnetic moment in the octupolar 
phase. Thick line: (M z ) staggered magnetization, thin line: 
(Tf) octupolar moment, as a function of the uniaxial pressure 
a || [100] (a in arbitrary units). 

Stress applied along the z-axis induces 0\ which trans- 
forms according to the identity representation Ai g , thus 
it does not appear in the invariants, and it is not pre- 
dicted to induce magnetism. 

We note that the A\ u triakontadipole J x J y J z (J x — J y ) 
(see Table [J would not give rise to stress-induced mag- 
netism and is therefore not a suitable choice as order 
parameter in the limit H — > 0. 

We now discuss the high-field behavior at H > H CIt ±. 
The single-ion levels t\ and d_ would cross at H CTOSS — 
37. 3T. Since \tx) and |d_) are connected by E opera- 
tors including O zx , a range of fields centered on H cmss 
is certain to favour {O zx ,O yz } quadrupolar order, and 
simultaneous {J x ,Jy} dipolar order. We chose a weak 
quadrupolar interaction A qua d = 0.054K in Eqn. J3J; 
this gives quadrupolar order between the critical fields 
H ctfl = 35.8T and H CIy3 = 38.8T. The amplitude of 
quadrupolar order is not small (Fig. ^) but the ordering 
temperature is low (~ IK) because the coupling is weak. 
The E phase shows up as the steep part of the magne- 
tization curve in Fig. 0] (right). For A qua( j = we would 
have a jump- like metamagnetic transition at H — H closs . 

In summary, our crystal field scheme gives an H-T 
phase diagram in overall agreement with experiments. 
Time reversal invariance breaking by the Tf octupolar 
order in the low-field phase is essential to allow the pre- 
diction of large-amplitude antiferromagnetism induced 



by transverse uniaxial stress. The disjoint high- field 
phase has mixed quadrupolar-dipolar character. 
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